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INTRODUCTION 


Dans la plupart des secteurs d'activités participant au processus de conception de structures, l'identification d'une 
structure optimale tient une place centrale. Cet optimisation porte généralement sur l'allègement des pieces tout 
en conservant certaines propriétés mécaniques telles que la rigidité ou la résistance du matériau. En pratique, deux 
contraintes mécaniques jouent un rôle déterminant pour la bonne tenue d'une structure. La première est relative à 
l'intégrité du matériau sous les efforts appliqués. La seconde contrainte mécanique révélatrice du comportement 
de la structure correspond à la résistance aux phénomènes dits de flambage. L'objectif de ce rapport est de 
présenter une algorithme d'optimisation topologique pour la conception de structures sous contrainte de flambage. 


Problématique et Solution : 
JUSTIFICATION DE CHOIX : 


Comment alors déterminer des structures innovantes avec la 


Dans le domaine de génie civil, l'identification minimisation de la masse afin de réduire les coûts de fabrication? 


de la meilleure structure selon divers critères Une approche possible est basée sur la modélisation mathématique 
mécaniques est prépondérante lors de du problème de conception en termes de problème d'optimisation et 
la conception de ponts, de bátiments,..... sur sa résolution efficace. En particulier, l'optimisation topologique 


permet de découvrir des dispositions de matières que l'humain 


n'aurait pas pressenti. 


Quelques 
définitions 


Au cours de l'optimisation , l'analyse du comportement de la structure peut être 
réalisée à travers une simulation numérique. L'une des principales techniques est 
basée sur la méthode des éléments finis, consistant à discrétiser la structure 
selon un type d'éléments qui peut être uni, bi ou tri-dimensionnel dont chaque 
extrémité est définie comme un noeud de la structure. Le résultat de cette 
discrétisation est appelée un maillage. Selon le type d'éléments finis utilisés, on 
distingue deux catégories de structures : d'une part les structures dites 

treillis, constituées d'éléments uni-dimensionnels, et d'autre part les structures 
dites continues, discrétisées à l'aide d'éléments bi et tri-dimensionnels. 





Le flambage ou flambement 
est un phénomène 
d'instabilité d'une structure 
élastique qui pour échapper à 
une charge importante 
exploite un mode de 
déformation non sollicité mais 
opposant moins de raideur à 
la charge. 


Optimisation 
topologique PRET 


-On distingue trois types d'optimisations de structures : l'optimisation dimensionnelle, l'optimisation de forme et 


l'optimisation topologique. 
-L'optimisation topologique consiste à rechercher, dans un espace de conception, la distribution optimale de 


matière représentant la structure. 

-On estime que, pour des problèmes de conception courants, l'optimisation d'un critère de performance peut 
mener à des gains variant de 40 a 100 % avec l'optimisation topologique.Tandis que, dans le cas de l'optimisation 
dimensionnelle, les gains oscillent entre 5 et 10 % et entre 10 et 30 % pour de l'optimisation géométrique . Ainsi, 
l'optimisation topologique prend une place prépondérante lors de la conception de structures, notamment parsa 
capacité à proposer des designs optimaux innovants. 

L'optimisation topologique a de nombreuses applications dans les domaines de l'aérospatiale, de la mécanique, 


de la biochimie et du génie civil. 





(b) Structure optimale 
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Analyse du comportement 
dé là structure 





Calcul dé nouvelles 
variables de conception 


| Nouvellestructure | 








Structure optimale * 





Optimisation topologique 
de structures sous 
contraintes de flambage 


nous portons notre attention sur l'optimisation de structures treillis, constituées de poutres, sous 
contrainte de flambage. L'objectif est de minimiser la masse en optimisant le nombre d'éléments 
structuraux et les dimensions des sections transversales associées. Dans un premier temps, nous 
explicitons l'analyse de flambage pour des structures à base de poutres. Ensuite, nous formulons le 
problème d'optimisation comme un problème d'optimisation non-linéaire en variables mixtes. Puis 
nous proposons une fonction objectif adaptée et un algorithme de résolution adapté au problème 
considéré. Enfin, sur un ensemble de problèmes issus de la littérature, nous illustrons les résultats 
numériques obtenus. 


Pour l'analyse du flambement Nous sommes donc en mesure d'étudier le problème aux valeurs propres généralisé lié au flambage 
correspondant à : 
(K+ AKo) = 0 
où les valeurs propres à correspondent à des facteurs scalaires de vecteur global des forces F, 
et les vecteurs propres v représentent les déformations de la structure correspondantes. La plus 
petite valeur propre positive Acr est appelée charge critique de flambage. Pour rappel, on dit 


qu'une structure résiste au flambage lorsque la valeur de Acr est supérieure a 1. 


Optimisation topologique 


de structures sous les éléments structuraux étant des poutres, une unique contrainte 
. est nécessaire pour s'assurer que la structure ne flambe pas 
contraintes de flambage globalement et localementsous les efforts appliqués : 
1 - Acr(x) < 0 


Cependant, avant d'observer un flambage local ou global, certaines zones de sections transversales présentant des épaisseurs fines 
peuvent flamber. Pour prévenir ce type de flambage, appelé aussi flambage de peau, nous définissons les contraintes pour un élément 
Ii: zi*hi(x)<0 
Avec: hi: une fonction dépend du type de la section transversale. 
zi: associée à chaque élément / pour indiquer sa présence dans la structure optimisée . Zi = 1 si i est présent . Sinon zi =0 


Minimiser f (x) 
x={z}€/0, LIN x(R+)Np 


1-Acr(x) < 0 :flambage globale 
zi"hi(x) <0 :flambage de peau 


zi € {0, 1}; =1,...,N 





la minimisation de la masse de la structure est souvent considérée comme 


Optimisation topologiq ue fonction objectif et est définie comme : 
de structures sous ON 
contraintes de flambage f(2)=2 pes A(d) 


nous proposons une fonction objectif où les propriétés géométriques de chaque élément, et pas 
seulement la masse, sont utilisées pour définir le critère d'optimisation. Plus spécifiquement, nous 
proposons de minimiser la masse de la structure ainsi que la somme des moments d'inertie selon les 
axes yet zZ. 





Minimiser f (x) 


Alors le probléme non linéaire de minimisation devient : 
Minimiser 
1 - Acr(x) <0 :flambage globale 
zi"hi(x) <0 :flambage de peau 
N 1 dmin < di < dmax 
fla): Y pr (1.A(d) y UoN Uuldi) + 1.(d;))) 


i=1 i=1 “1 


zi € {0, 1}; =1,...,N 





Algorithme de résolution 


A l'état initial, tous les éléments sont considérés présents dans la structure (cela correspond à z = 1). Dans un premier temps, les valeurs 
des variables z sont fixées et le problème avec les variables continues d est résolu. C'est un problème non-linéaire . donnant les dimensions 
d optimales des sections transversales des éléments présents dans la structure. Dans une seconde étape, à partir de la solution .nous 
traitons de nouveau les variables binaires. Estimer les valeurs des variables binaires z revient à étudier la suppression des éléments 
structuraux (fixer zi = O correspond à la suppression de l'élément i). Ici, une méthode pour supprimer les éléments un par un est présentée. 
Lorsqu'un élément structural est enlevé, une nouvelle disposition de poutres est obtenue et les deux étapes précédentes sont répétées. 
Afin de déterminer la vrai influence de chaque élément, nous devons déterminer les dimensions optimales des sections transversales 
associées. Pour rappel, la résolution du problème non-linéaire pour une configuration de structure donnée, nous donne les dimensions d 
optimales des sections transversales. À partir de ce résultat, un ensemble S de poutres, dont les aires des sections transversales sont plus 
petites ou égales à un seuil AT, peut être défin : S := {i | A(di) < AT } 

Ensuite, pour chaque élément i de S, nous estimons l'impact de la présence de l'élément ipar une analyse de sensibilité sur la charge critique 
de flambage Acr: ai:= AÀcr(d, z1, ..., zi = 0, ..., zm) - Aer(d, z1, ..., zi = 1, ..., zm), Vi E S 

où le premier et le second terme correspondent respectivement à la charge critique de flambage de la structure lorsque l'élément i est 
absent et présent. Alors l'élément e, qui maximise ai, est supprimé de la structure (i.e. nous fixons ze = 0) Avec : e := argmax ai 

Une fois qu'un élément a été supprimé, nous obtenons une nouvelle structure sur laquelle on répète le même procédé . A l'état initial, la 
valeur du seuil AT est égale à l'aire minimale définie par les bornes sur les dimensions des sections transversales. 


Algorithme de résolution 





La procédure proposée est présentée dans l'Algorithme : (P”) , 

input: dmin, dmax, o etp Minimiser f (x) 
Initialisation : z = {1}m et AT. 
Résoudre (P”) 


Calcul de la masse structurale m selon z et d 
Mx*=mM;Z*=Z et d*-— d 


1 - Acr(x) <0 :flambage globale 
for j = 1,2 do zi"hi(x) <0 :flambage de peau 
Construire S := {i | A(di) <j x AT] dmin < di < dmax 
while S +10 do 


Calcul des sensibilités ai - i B 
Déterminer e = argmax aii € S zi € {0, 1); v1=1,....N 
Fixer ze — 0 
Résoudre (P”) 
Calcul de la masse structurale m selon z et d 
if m < m=» then 
Mx=M;Zx*=Z et d*-— d 
end if 
Construire S := {i | A(di) <j x AT } 
end while 
end for 
output: mx, zx et dx 






Masse structurale optimale (kg) et temps de calcul (s) pour 


A | C at ONS à d ive rS chaque problème et chaque section transversale 


Section transversale creuse 
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TABLE 4.9 — Nombre de variables continues el entières à l'état initial 
Variables continues selon la section transversale 
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(a) Structure initiale 
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(c) Mode de flambage [d] Evolution de la masse strnéturalé pour les diferentes sections transversales 


l'optimisation topologique est utilisable aussi pour les structures continue a | aide Méthode 
du gradient topologique : la notion de gradient topologique quantifie l'influence sur un 
critère d'optimisation de l'insertion d'un trou dans une structure l'état initial, le 
domaine de conception est considéré rempli de matière .On dit que la fraction de 
volume a occupé est de 1. Ensuite, l'objectif est de supprimer successivement des 
fractions de volume Aa. à partir des champs des gradients topologiques. 


Elargissement 






Structure initiale O 
Définition de N, a, Ag, Amin et AG; 





Calcul de T pour la 
structure ( 


Gala): 







T tel que [N°] = [N (æ — Ax)] 


16 (9) — Gs (Q°)| 


Ia <=“ 


-AGO est le seuil d'acceptation 
de variation de la fonction Go. 


-Le domaine Or correspond aux 
éléments dont la valeur du 
gradient topologique est 
supérieure . 


N 
l ] 3 
In E apra) <0, où plej = (22 - 1) 
N T 
4 pan | 













(a) Structure N (b) Structure Ng 


FIGURE 5.4 — Concept du gradient topologique 





(a) Champ T des gradients topologiques engendré (b) Champ T des gradients topologiques engendré 
05 par = 0.1 


par Aer = O. 
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(d) Domaine engendré par Aa = 0.1 


(c) Domaine engendré par Aa = 0.05 


Logiciels de l'optimisation Topologlque 





L'optimisation topologique permet de déterminer la forme et la répartition de matière la plus efficace 
en fonction du matériau défini et des contraintes. La forme générée, tout à la fois plus légère et plus 
rigide, a un aspect organique et nécessite moins de matière. Ceci est un enjeu majeur dans les 
domaines de l'aéronautique, le médical, l'automobile,Genie civil ou la machine spéciale. Generative 
design est la réponse dont les bureaux d'études, les designers et les architectes ont besoin afin de 
concevoir rapidement et facilement des produits efficaces grâce à une nouvelle façon de concevoir 
des formes. 


Exemples des logiciels de | OT 
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Application 
au Maroc 








-Dans le Maroc, comme dans le monde entier, l'optimisation topologique des structures est une 
technique fréquemment utilisé et qui a 

trop évolué dans ces dernières années Parce que de nombreux logiciels ont été développés dans 
ce sens. 

-Cette technique est fortement utilisée dans le secteur de Genie civil et de mécanique au Maroc, 
beaucoup de bureau d'études de fabrication utilisent l'optimisation topologique à l'aide de 
solidworks et Robobat afin d'avoir une structure performante avec un minimum de masse . 


Conclusion 











Dans cette présentation, nous avons étudié un problème d'optimisation 
d'une structure, sous des contraintes de flambage. Nous avons présenté : 
— À travers une formulation non-linéaire en variables mixtes, nous avons 
proposé une fonction coût reposant sur une somme pondérée de la masse et 
des moments d'inertie de chaque élément. 

— Un algorithme a été développé, basé d'un côté, sur un problème non- 
linéaire pour optimiser les dimensions des sections transversales des 
poutres, et de l'autre sur un critère utilisant une analyse de sensibilité de la 
charge critique de flambage pour supprimer un par un les éléments de la 
structure. 

-Les résultats numériques, sur un ensemble d'exemples de la littérature, 
montrent que l'augmentation du nombre de degrés de liberté décrivant l'aire 
des sections transversales mène à des réductions de masse structurale 
jusqu'à 7,15%. 
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